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discussed. 
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I. INTRODUCTION 



The existence of observable effects originating from the quantum nature of the electromag- 
netic field has received much attention since 1948 when Casimir predicted that zero-point 
field fluctuations give rise to an attractive force between two neutral conducting plates at 
rest in the vacuum [1]. The same year Casimir and Polder provided an explanation for the 
retarded long-range van der Waals interaction between two neutral polarizable objects as a 
manifestation of the zero-point energy of the electromagnetic field [2]. They found that re- 
tardation yields a decay law of the interaction energy as r -7 at large interatomic separation 
(Casimir- Polder potential). Casimir's results also showed how geometrical constraints can 
affect vacuum field fluctuations. These effects have been measured experimentally, and the 
results obtained are in good agreement with the theory [3-5]. More recently the attention of 
theoreticians has been also drawn by how a dynamical change of geometrical or topological 
boundaries affects vacuum field fluctuations, giving rise to observable effects such as modi- 
fications of the Casimir effect [6] or creation of real quanta from the vacuum (the so-called 
dynamical Casimir effect) [7, 8]. Dynamical Casimir effect is closely related to the Unruh 
effect, which establishes that an atom or a charge uniformly accelerated in the vacuum be- 
haves as if it were immersed in a bath of thermal radiation with a temperature proportional 
to its acceleration [9]. The concept underlying all these phenomena is that the notion of 
vacuum and its physical properties depend critically on the physical system considered and 
on the boundary conditions. 

Initially, the interatomic Casimir-Polder (CP) potential was understood in terms of the 
energy of the zero-point fluctuations of the electromagnetic field. More recently, it was shown 
that it can be also obtained as a consequence of the existence of correlations between the 
fluctuating dipole moments of the atoms, induced by the spatially correlated vacuum fluctu- 
ations. In other words, the CP interaction energy between two atoms in their ground state 
can be seen as the classical interaction energy between the instantaneous atomic dipoles, 
induced and correlated by the spatially-correlated vacuum field fluctuations [10, 11]. This 
model is conceptually intriguing because gives a classical picture of CP forces: the quantum 
nature of the electromagnetic field enters only in the assumption of vacuum fluctuations 
as a "real" field affecting atomic dynamics. This model has been also generalized to the 
three-body CP potential between three atoms in their ground or excited states [12, 13]. In 
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this case, any pair of atoms interacts via their dipole moments which are induced and corre- 
lated by the vacuum field fluctuations, modified (dressed) by the presence of the third atom. 
Because the presence of one atom modifies the spatial correlations of the electric field, the 
interaction between two atoms changes if a third atom is present, and this evantually yields 
a non-additive interaction. Thus CP forces between atoms are a direct manifestation of the 
existence of nonlocal correlation of zero-point fluctuations, and so their measure can be used 
as an indirect evidence of field correlations and for investigating their nonlocal properties. 

Many conceptual difficulties in quantum mechanics are involved in the notion of nonlocal 
correlations, also in connection with relativistic causality. In particular, Hegerfeldt reopened 
the question of nonlocality and causality in the QED context on a quite general basis [14- 
20]. Specific calculations have shown that the dynamics of local atomic or field operators 
is causal, but that the correlation of atomic excitations of two spatially separated atoms 
exhibits a nonlocal behaviour, being different from zero even if the two atomic sites have 
a spacelike separation [17, 21]. Non-local terms appear also in the spatial correlations of 
the energy density during the dynamical dressing/undressing of a static source interacting 
with the relativistic scalar field [22]. The question of relativistic causality and its relation 
with the nonlocal correlations of vacuum fluctuations has been also examined in connection 
with the Unruh effect [23]. All this emphasizes the conceptual importance of investigating 
whether nonlocal correlations of vacuum fluctuations may be at the origin of observable 
effects. Recently we have examined the question of causality in the dynamical CP interaction 
between an excited and a ground-state atom during the dynamical evolution of the excited 
atom [24]. If one of the two atoms is in the excited state at time t — 0, the interaction 
energy between the two atoms is non- vanishing only after the causality time t = R/c, R 
being the interatomic distance. This indicates that causality is a basic property of QED as 
far as local quantities such as field energy densities or the two-body CP forces (which are in 
principle measurable by a local observation) are considered. 

A question worth considering is what happens when nonlocal quantities, such as the 
three-body forces, are considered. This is related to the very nature of many-body forces, 
which appear to be inherently nonlocal because they cannot be measured through a single 
local measurement. We have recently investigated this issue in the time-dependent three- 
body CP interaction between three atoms initially in their bare ground-state, during their 
dynamical self-dressing [25]. Our results indeed indicate that there exist time intervals and 
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geometrical configurations of the three atoms for which the three-body interaction energy 
exhibits a nonlocal behaviour, related to nonlocal properties of the field correlation functions. 
This should allow to investigate the nonlocal properties of the field by measurements of the 
(observable) three-body Casimir-Polder forces. 

In this paper we address a similar question for the dynamical CP potential between one 
excited- and two ground-state atoms, during the dynamical self-dressing of the excited atom 
and the initial part of its spontaneous decay. We investigate the causality problem in the 
time-dependent three-body CP potential and its relation to the nonlocal properties of the 
field emitted by the excited atom during its short-time evolution. Compared to the case of 
three ground-state atoms, a resonant contribution is now present both in the field emitted 
by the excited atom and in the time-dependent potential; this new contribution arises from 
an additional term in the correlation of the induced dipoles of the two ground-state atoms. 

The paper is organized as follows. In section II we consider a two-level atom (say C) 
initially in its excited state and we investigate the dynamics of the spatial correlations of 
the electric field during its dynamical evolution. We show that this correlation has a non- 
local behaviour. In Section III we calculate the interaction energy between a pair of atoms 
located at some distance from atom C and show that this energy has nonlocal properties as 
a consequence of field nonlocality. Finally, in section IV, we calculate the total three-body 
CP potential by an appropriate symmetrization of the role of the three atoms and discuss 
how the nonlocal behaviour of the field correlation function influences the dynamical CP 
potential between the three atoms. 

II. THE FIELD CORRELATION FUNCTION 

We first evaluate the correlation function of the electromagnetic field during the sponta- 
neous decay at short times of a two- level atom (C), initially in its bare excited state. We 
describe our system using the multipolar coupling Hamiltonian, in the Coulomb gauge and 
within the dipole approximation [26], 

H = H C + H F + H mt (1) 

where 
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where o;o is the atomic transition frequency, S z , S± are the pseudospin operators of atom C, 
the bosonic annihilation and creation field operators and Yc is the position of 
atom C. Ckj is the coupling constant given by 



where /x c is the matrix element of the electric dipole moment of atom C. 

The initial state is assumed as the factorized state 
| vac, with the atom C in its bare excited state and the field in the vacuum state. First 
we wish to evaluate on the initial state | vac, |c) the average value of the equal-time spatial 
correlation of the field at two different points and r#, (d_u(rA,t)d± m (rB,t)) (we work in 
the Heisenberg representation), where 



is the transverse displacement field operator (the momentum conjugate to the vector poten- 
tial, in the multipolar coupling scheme) which, outside the atoms, coincides with the total 
(transverse plus longitudinal) electric field operator [27]. From now on we shall use the 
symbol E in place of dj_. 

Our approach closely follows that used by Power and Thirunamachandran in [28, 29]. 
Solving Heisenberg equations of motion for the field operators a^j(t) and a^^i) at the second 
order in the coupling constant, we obtain the following expansion for the field operator 
[28, 29] 



The operator E^(r, t) is the free-field operator at time t, while E^(r,t) and E( 2 )(r, t) 
are source- dependent contributions. Explicit evaluation of (5) shows that both EW(r,t) and 
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(3) 





E(r, t) = E (0) (r, t) + E (1) (r, t) + E (2) (r, t) 



(5) 



E( 2 )(r, t), contrarily to E^(r, t), contain the Heaviside function 6(ct—R), where R —\ r—r c 
is the distance of the observation point r from atom C, 



This expresses a causal behaviour of the source electromagnetic field. Hence the electric 
field at the second order can be expressed as the sum of two terms: a free-field contribution, 
which is independent of the presence of atom C, and a source-dependent contribution which 
is strictly causal, 



It should be stressed that these results have been obtained in the multipolar coupling 
scheme, where the operator conjugate to the vector potential is the transverse displacement 
field; outside the atoms, it coincides with the total electric field, which obeys a fully retarded 
wave equation. In the minimal coupling scheme, on the contrary, the conjugate momentum 
is the transverse electric field, which obeys a wave equation with the transverse current 
density as source term; in this scheme we would have obtained a non-retarded solution, and 
electrostatic terms should be added in order to restore a causal propagation of the field. This 
illustrates the remarkable advantage of using the multipolar coupling Hamiltonian, which 
is obtained from the minimal coupling Hamiltonian by the application of the Power-Zienau 
transformation [27]. We now evaluate the expectation value of the correlation function of 
the electromagnetic field {vac, jcl Ee(r A ,t)E m (r B ,t) \ v ac, fc) & t the two points r A and 
r£. Up to the second order in the electric charge e and using the expressions for the field 
operator given in [28, 29], this correlation function is obtained as 



E (1) (r,t) ~6(ct-R); E (2) (r,t) ~ 6(ct - R) 



(6) 



E(r,t) = E (/ree) (r,t) + E ( ~°(r,t) 



(7) 



(E e (r A ,t)E m (r B ,t)) = (E^(r A , t)E%\r B , t) 
+E { l\v A , t)E£>(r B , t) + Ef \r A , t)E%\r B , t) 
+E$>\r A ,t)E£>(r B ,t)) 



(8) 



with 



(Ef\r A ,t)E£\r B ,t)) 



2ixhc 



£(^)i(eitf)mfce*-^- r *>; 



(9) 



V 
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{Ef\v A} t)E^(v B ,t)) = tftfFi e —F« p e -—6(ct - 0)9{ct - a); (10) 

(E?\r A ,t)EV>(r B ,t) + Ef\r A ,t)E%\r B ,t)) = ^F^-F^-^-9{ct - P)9(ct - a) 
~y {E ^nV?(e%)Ke k ,)ne- ik - R -^F« p i(e-^ - e *« e -<*+ fc >*)6>(ct - 0) 

+ (c.c.(A# B,a^(3,(ln) # (mp)))| (11) 

where a =| r B — r c | and /3 =| — r c \- 

Eq. (9) describes the zero-point contribution to the field correlation function and does 
not play any role in the causality problem for three-body forces we are concerned with. On 
the other hand, the two terms (10) and (11) depend explicitly on the position of atom C. In 
particular, the term (10) arises from the retarded field emitted by atom C at the two points 
i% r_B and it is causal. This is expected since the electric-field operator d^(r, t) vanish for 
t < r/c. The other contribution, which contains both the second-order field (causal) and 
the free-field at time t, is responsible for the non-local behaviour of the field correlation 
function. In order to discuss this point, we partition the field correlation function in the 
following form (disregarding the free-field contribution which, as mentioned, is not relevant 
for our purposes) 

(vac, Tc| E £ (r A ,t)E m (r B ,t) | vac,] c ) 
= (E e (r A ,t)E m (r B ,t)) nr 

+(E e (r A ,t)E m (r B ,t)) r (12) 



where 



{E e {r A ,t)E 

in ( r B,t)) nr 

= -y{E^Vf(e k ,)Ke kj )ne- k - R - 

I kj 

x VTk F ™ p ^ e ~ ika ~ eikoae ~ i{ho+k)ct ) d ( ct - P) 

+ (c.c.(A^B,a^ P,(ln) # (mp)))\ (13) 



is the nonresonant contribution to the correlation function, and 

(E e (r A ,t)E m (r B ,t)) r 



- * tripsin ^ ^ _ PJ 

x0(ct-a) (14) 

is the resonant contribution, which derives from the pole at k = ko in the frequency inte- 
gration. The nonresonant term (13) is equal but opposite in sign to that already obtained 
when atom C is in the ground state [25]. The resonant term is not present in the case of a 
ground state atom, of course. Inspection of (13) and (14) clearly shows that if the two points 
ta and r B are outside the causality sphere of atom C, that is if a,/3 > ct, the correlation 
function (12) reduces to zero. When both points and r# are inside the light-cone of atom 
C, the correlation function is modified by the presence of atom C. All this is compatible 
with relativistic causality, of course. Yet, nontrivial results are obtained if just one of the 
two points and r B is inside the causality sphere of atom C. For example, when a < ct 
and f3 > ct the correlation function is modified by the presence of atom C. Moreover, this 
happens whatever the distance between the two points and v B . This result indicates non- 
local features of the field correlation function, which originate only from the non-resonant 
part of the correlation function, as clearly shown by Eqs. (13-14) 

III. THE THREE-BODY CONTRIBUTION TO THE DYNAMICAL CASIMIR- 
POLDER INTERACTION 

Let us now consider two more ground-state atoms, A and B, located at points and r# 
respectively. We wish to evaluate their Casimir-Polder interaction energy AE C (A, B) in the 
presence of atom C. Our aim is to investigate whether the nonlocal behaviour of the field 
correlation function discussed in the previous Section may reveal itself in the time- dependent 
interaction energy between the two atoms. This is indeed expected because it is known that 
the Casimir-Polder interaction between two atoms depends on the vacuum field correlations 
evaluated at the atomic positions [10]. We have already discussed a similar problem in the 
case of three atoms initially in their bare ground state [25]; the main difference in the present 
case is the presence of a resonant contribution to the correlation function. Our approach is a 
generalization to the time-dependent case of the model already used to calculate the three- 
body potential with one atom excited in a time-independent approach [13]. Following the 
same arguments used in [13], to which we refer for more details, the three-body contribution 
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to the interaction energy between atoms A and B in the presence of the excited atom C 
consists of two terms. The first is related to the non resonant part of the correlation function 
and is formally equivalent to that obtained when the three atoms initially are in their bare 
ground state. The second is related to the resonant part of the field correlation function. 
Thus we write the interaction energy between A and B in the presence of C as 

AE C (A, B) = AE C (A, B) nr + AE C (A, B) r (15) 

where the first term is a non-resonant contribution and the second the resonant one. These 
two contributions are expressed as 

AE c (A,B) nr = ]T a A (k)a B (k') 
x (E e (kj,r A ,t)E m (k'j',r B ,t)) n . r y lm (k,k', 1 ) (16) 

and 

AE C (A, B) r = ]T a A (k)a B (k') 
kjk'i' 

x (E e (kj, r A , t)E m (k'j f , r B , t)) r V[ m (k , 7). (17) 

Ee(kj, r, i) are the Fourier components of Et(r,t), 

V lm (k,k',r) = —F^-icosk^ + cosk'j) (18) 
2 7 

is the classical potential tensor between oscillating dipoles at frequencies k and k! [12], and 

VL{h,i) = -F2 m (^^ (19) 

is the potential tensor for dipoles oscillating at the resonant frequency k . 7 =| — r B \ 
is the distance between dipoles A and B, and Fj n = (—V 2 5e n + V^V„) 7 is a differential 
operator acting on the variable 7. The resonant contribution(17) is specific to the excited- 
atom case and does not appear when all atoms are in their ground state. After lengthy 
algebraic calculations, we obtain the explicit expressions of AEc(A, B) nr and AEc(A, B) r 



AE C (A, B) nr = ^^F^FlF^-L- f {(^ k(J3 + 7) + sin k(J3 - 7)) 

xe- lka 9(ct -a) + sin k{3 [ e - ifc(a+7) #(ct - (a + 7)) + sgn(a - -f)e- tkla -^9(ct- \ a - 7 |)" 
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+c.c. (a # 

+^^ 2 /'J 2 ^»Jl^^;{ai>(*o)e- i ™ , [e-*°e(ct - a) 

Jo A; + A; 

+ ( e - ife ° (a+7) 0(c* - (a + 7)) + sgn(a - ^e'^^eict- \ a - 7 |)) J™ dka A (k)^^e- iuJkt 

+a A (k )e i ** t (c.c.(a^P))\ (20) 

AE C (A, BY = W« A (A;o)aB(A;o)2K(F/ n F- p ^^)F7 m (-^^)^(ct - ff)0(ct - a) (21) 

In order to investigate possible evidence of nonlocality in the three-body Casimir-Polder 
interaction (15), let us consider a few limiting cases. For a <C ct and (3 <C ct and in the limit 
of large times (compatibly with the perturbative expansion we have used), the interaction 
energy between A and B in their ground states reduces to the value 

AE C (A, B) = -^tfF? m Ff n FZ p — 



dk a A (k)a B (k) 



a/?7 

cos ka( smk(3 + 7) 
k + k L V 

+ sin k{(3 — 7)) + sin k(3 ( cos k(a + 7) 
+sgn(a — 7) cos k(a — 7)) + (c.c.(a; # (3^j j 

+tftfa A {k )a B (k )F°^ n Fj m ^ 

x (cosA;o(a; - (3 + 7) + cos fc (a — — 7)) (22) 

which is already known from time-independent calculations [13]. This means that, after a 
certain time, the interaction energy settles to a quasi- stationary value, as expected [24]. 

A noteworthy result is obtained when we consider a time t such that a > ct and/or (3 > ct. 
This means that at least one of the two atoms A and B is outside the causality sphere of 
C. Quite unexpectedly, equations (20) and(21) show that in this case the interaction energy 
between A and B is affected by the presence of C. In order to point out the most relevant 
aspects, let us focus on the specific configuration a, (3 > ct and 7 < ct, that is A and B 
outside of the light cone of C but inside the light cone of each other. This configuration of 
the atoms and their causality spheres are schematically illustrated in Fig. 1. The interaction 
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energy AE C (A,B) 



is then 




Jo k + k 



- - e -iko(P-i) _ _ sgn (/3 - 7 - ct)e ikoW - y) 
+ (c.c.(A^B,a^/3) )\ 



(23) 



The main point is that there are time intervals for which the expression above does not 
vanish: this happens when ct < a < 7 + ct and/or ct < j3 < 7 + ct. We stress that in such 
cases both atoms A and B are outside the light come of C: nonetheless their Casimir-Polder 
interaction energy is affected by atom C, indicating nonlocal aspects in their interaction 
energy. This does not contradict the fact that in this case the correlation function (12) 
can be zero if both A and B are outside the light-cone of C. In fact, the calculation of 
the quantity AE C (A,B) involves a sum over the field modes of a product of the electric 
field Fourier components and of the interaction potential Vi m , which also depends on k. We 
also observe that this effect derives exclusively from the non-resonant contributions to the 
three-body CP potential: the resonant three-body CP potential is non-vanishing only when 
both atoms A and B are inside the causality sphere of C. Thus it seems that the nonlocal 
properties of the electromagnetic field emitted by atom C during its dynamical self-dressing 
become manifest in the time-dependence of the three-body CP interaction energy between 
atoms A and B, and only the (nonresonant) virtual processes contribute to this effect. 

An important conceptual point is the physical meaning of the interaction energy 
AE C (A, B). It is not a potential energy related to a single atom or to the whole sys- 
tem of the three atoms, but it is related to the change of the interaction between two atoms 
(A and B) due to the third atom (C). Therefore its measurement must necessarily involve 
some correlated measurements on both atoms A and B, in order to separate it from other 
contributions to the three-body energy such as AE B (A, C) and AEa(B, C). 
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FIG. 1: A configuration of the three atoms A, B, C at time t such that a > ct, (3 > ct and 
7 < ct, for which a nonlocal behaviour of the interaction energy AEq(A,B) is found. Red, blue 
and green (long-dash-dot, dashed and continuous lines, respectively) circumference's arcs of radius 
ct represent respectively the causality spheres of atoms A, B and C at time t. 

IV. THE TIME-DEPENDENT THREE-BODY CP POTENTIAL BETWEEN 
ATOMS A,B AND C 

We now evaluate the following quantity, obtained by a symmetrization of the interaction 
energies of any pairs of atoms in the presence of the third one 



AE(A, B, C) 




x (EtQtj, i\4, t)E m (k'j', r B , t)) nr V Em (k, k', 7) 
+(A -> B -> C) J 



+ (yA(k)a B (k')(E i (kj,r A ,t)E m (k'j',r B ,t)) r 



kjk'j' 



xVf m (k ,j) 



(24) 
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where (A — > B — > C) indicates terms obtained from the first double sum by a permutation of 
the atomic indices. In stationary cases this quantity has been shown to be equivalent to the 
three-body Casimir- Polder potential, as obtained by sixth-order perturbation theory [13]. 
This motivates our choice to consider this physical quantity. We stress that we symmetrize 
only on the nonresonant part, for which the role of the three atoms is indeed symmetrical; 
the resonant part should not be symmetrized because the contribution of the three atoms 
to this term is not symmetrical, only C being in an excited state. 

We now wish to investigate if the nonlocal aspect discussed above for the interaction 
energy AE C (A, B) are present in AE(A, B,C) too. Explicit evaluation of (24), yields 



AE(A, B, C) = AE {I) (A, B, C) + AE {II) (A, B, C) 
+AE {r) (A,B,C) 



(25) 



where 



he f°° \ / 

AE (I) (A, B, C) = —F^F^— 1 1 dua A (iu)a B (iu)ac(iu) \e~^^ [6 

— sgn(a — ct) — sgn(/3 — ct) — sgn(7 — ct) — sgn(o; + (3 — ct) — sgn(o; + 7 — ct) 

-sgn(7 + (5 - ct)) 

+ e -«(«-0+r) (- S gn((3 - ct) + sgn(a + 7 - ct)) + e- u{a+ ^ (-sgn( 7 - ct) 
+sgn(o; + (3 — ct) 



+ e - u{ ~ a+l3+l) (-sgn(a - ct) + sgn(/5 + 7 - ct) 



(26) 



and 



he if 

AE (II) (A,B,C) = —F&F^— ^a B (k ) 



cosko(a-ct) / duaA(iu)ac(iu) 
Jo 



-u\j3— 7— ct\ 



x (V M(/3+7+ct) + sgn (/3 + 7 _ ct ) e -«|/3+ 7 -rf| + sgn^ _ 7 _ ct y 
+sgn(/? - 7 + ct)e- M|/3 ~ 7 + rf| ) 

' k ° {a " Ct) / o °° du u a A {iu)a c {iu) ( e -P+T**) + e -«l/»^--l + e -«l/>-7-*| _ e -«l/»-rK*| 

cos fc (/3 + 7 - ct) a B (m)a c (iu) (e- M(Q+rf) 



+ ( S # C, /? # 7 ) 



9(ct - a) + 4a A (/c ) 



+ sgn(a - ct)c- u \ a - ct \) + sinfc o(/^ + 7 ct) duu aB{m)ac{m) L-u (a+ ct) + e - u |a-ct| 
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x0(ct-(/3 + 7 )) 

+ (A^B, a^(3) + (A^C, a# 7 )|; (27) 
are the nonresonant contributions, while 

A£(A, 5, C) (r ) = -/^fMMMW^(F/ n F« p — -^)0( c * - 0)0(ct - a) 

(28) 

is the resonant one. We have assumed isotropic atoms, that is = ||^| 2 ) and ot{iu) is 
the dynamic polarizability extended to imaginary frequencies. 

Eq. (25) describes the time-dependent symmetrized three-body CP potential as a function 
of time for a generic configuration of the three atoms (for times shorter than the spontaneous 
decay time of the excited atom, due to the limitations of our perturbative treatment). As 
in the case discussed in the previous Section, we now consider specific cases relevant for 
our discussion about nonlocal aspects of the dynamical interaction energy. First of all, it 
is immediate to see that AE(A, B, C) vanishes if each atom is outside the light cone of 
the other two, that is for a, (3, 7 > ct. On the contrary, AE(A, B,C) is non-vanishing 
for times such that each atom is separated by a time-like interval from the other two. In 
particular, for large times, the time-dependent terms rapidly decrease to zero and we find the 
well-known stationary result [13]. This means that after a transient characterized by a time- 
dependent interaction, the three-body interaction energy settles to the time-independent 
Casimir-Polder interaction between three atoms with an excited atoms. All these results 
are compatible with relativistic causality and similar to those previously obtained for the 
dynamical three-body Casimir-Polder interaction between three ground-state atoms [25]. 
However, when the spatial configuration of the three atoms is such that two of them are 
separated by a time-like distance we find a non- vanishing three-body interaction, even if the 
third atom is outside their causality sphere. For example, when the separation of A from 
the other two atoms is space-like, Eq.(25) yields 



AE(A, B, C) = -Fl v FiF> m — jf dux A {iu)a B {iu)a c {iu){e-< a+ ^ + e -«(-«+^)) 
-^n^KK F ^^{ *c(ko)e- ik °^ f dka A {k f U l^\ ^ (29) 
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which in general does not vanish. Thus the nonlocal features of field emitted by the atoms 
during their self-dressing are evident also in the three-body interaction energy AE(A, B, C), 
with features which may differ from to those of AE C (A, B). Eq.(29) shows also that the 
nonlocal features of AE(A, B, C) stem from the non- resonant contributions, so that they 
are exclusively due to the virtual photons dressing the atoms, as we have recently discussed 
in the case of ground-state atoms [25]. 

V. CONCLUSIONS 

We have considered the Casimir- Polder interaction energy between three atoms with one 
atom initially in its excited state, using a time-dependent approach. We have discussed 
the problem of relativistic causality in the interaction between the atoms and its connection 
with the non-locality of spatial field correlations. The spatial correlation function of the field 
emitted during the spontaneous decay of the excited atom has been first obtained. We have 
shown that a non-local behaviour appears, in agreement with previous results, and that it is 
related to a non-resonant contribution related to the emission of virtual photons. We have 
shown that a non-local behaviour appears also in the dynamical Casimir-Polder interaction 
between two other ground-state atoms, during the initial stage of the spontaneous decay of 
the first atom. We have suggested that the appearance of this non-local behaviour can be 
ascribed to the non-locality of the field correlation function and that this new phenomenon 
should be observable. Thus we conclude that the nonlocal properties of the electromagnetic 
field emitted by the atoms during their dynamical self-dressing may become manifest in 
the time-dependence of the Casimir-Polder potential. We remark that previous studies 
of causality in the time-dependent two-body Casimir-Polder interaction have not shown 
indications of non local behaviour [24]. Hence the causality problem appears quite more 
complicated and subtle in the case of the time-dependent three-body Casimir-Polder energy, 
where nonlocal aspects may become manifest. 
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